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ON A SITE OF X-RAY EMISSION IN AE AQUARII 
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ABSTRACT 

An analysis of recently reported results of XMM-Newton observations of AE Aqr within a hypothesis 
that the detected X-ray source is located inside the Roche lobe of the white dwarf is presented. I show 
this hypothesis to be inconsistent with the currently adopted model of mass-transfer in the system. 
Possible solutions of this problem are briefly discussed. 

Subject headings: accretion, accretion disks - X-rays: binaries - (stars:) white dwarfs - stars: 
individual (AE Aqr) 



1. INTRODUCTION 

AE Aqr is a close binary system with an orbital period 

o o 

9.88 h, mass ratio 0.58-0.89, and an inclination i ~ 55 ±7 
(Welsh, Horne, & Gomer 1995). The distance to the 
system is about 100 pc (Friedjung 1997). The degener- 
ate companion is a magnetized white dwarf rotating with 
a period of Pg — 33 P33 s and braking with a rate of 
Pq — 5.64 X 10~^^ss~^, which implies the spin-down power 
of (de Jager et al. 1994; Welsh 1999), 

L,d = 6 X 10^3 p-3 1,0 {P/Po) ergs-\ (1) 
Here /50 is the moment of inertia of the white dwarf ex- 
pressed in units of 10^°gcm^. 

The normal companion is a K3-K5 red dwarf which 
overflows its Roche lobe and loses material through the 
LI point towards the white dwarf. This material man- 
ifests itself in a form of the optical/UV continuum and 
emission lines. It is neither accreted onto the surface of 
the white dwarf nor stored in a disk around its magneto- 
sphere. Instead, it interacts with the magnetic field of the 
white dwarf via a drag term (King 1993) and is leaving 
the system without forming a disk (Wynn, King, & Horne 
1997; Welsh, Horne, & Gomer 1998; Ikhsanov, Neustroev, 
& Beskrovnaya 2004). 

The system X-ray emission has recently been stud- 
ied with XMM-Newton (RGS, EPIC) by Itoh, Ishida, 
& Kunieda (2005). As they have reported, the num- 
ber density of plasma responsible for the detected X-rays 



the white dwarf. According to their scenario the heat- 
ing occurs at an adiabatic shock located at a distance of 
^ 10^*^ cm from the white dwarf and the observed emission 
is powered by the gravitational energy of the streaming 
material. 

In this letter I show that the above mentioned hypothe- 
sis is inconsistent with the currently adopted model of the 
mass transfer in the system and thus, a question about a 
location of the source of X-ray emission in AE Aqr remains 
open. 

2. CAN THE X-RAY SOURCE BE ASSOCIATED WITH AN 
EXPANDED STREAM ? 

Within the hypothesis suggested by Itoh et al. (2005) 
the flow of hot material responsible for the observed X- 
rays is feeded by the stream flowing into the Roche lobe 
of the white dwarf through the LI point. The numerical 
simulations of Ha Doppler tomogram of AE Aqr (Wynn 
et al. 1997; Welsh et al. 1998; Ikhsanov et al. 2004) have 
shown that a distance to which the stream could approach 
the white dwarf, r^in, significantly exceeds its corotation 
1.5 X 10^ Mq^^ Pt!^ cm. This distance, in 



radius Tcor — ^.u ^ -lu -""q.s ^ 33 
the general case, is limited to rm 
TA is the Alfven radius of the white dwarf, and 



ro ~ 10 
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> max[ro,rA], where 



(2) 



1.8 X 10" cm 



~ (2 — 3) X 10 " cm. The observed spectrum has 
been well fitted using the 4-temperature (0.14, 0.59, 1.4 
and 4.6 keV) VMEKAL model (Mewe, Kaastra & Liedahl 
1995). The analysis of the centroids of N and OLya lines 
shows no evidence for any significant orbital Doppler mod- 
ulation. The widths of these lines are close to 1000 km s^^. 

Analyzing these results Itoh et al. (2005) have discarded 
a possibility that the detected X-rays are emitted from the 
surface of the white dwarf indicating that is a few orders 
of magnitude smaller than corresponding conventional es- 
timates in the post-shock accretion column and exceeds 
the radius of the white dwarf by almost two orders of mag- 
nitude. Instead, they have suggested a hypothesis in which 
the detected emission is associated with a heating and ex- 
pansion of the material streaming through Roche lobe of 
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is a distance to which the material could approach the 
white dwarf if its angular momentum along the whole tra- 
jectory remains constant (see e.g. Eq. 2.14 in Warner 
1995). Here q and a are the mass ratio and orbital sepa- 
ration of the system components. 

The radial distance represents a minimum possible 
distance to which the material flowing through the LI 
point can approach the white dwarf. The value of this 
parameter does not depend on either the magnetic fleld 
strength of the white dwarf or the structure of the inflow- 
ing material and is only based on the angular momen- 
tum conservation law. For the material to come closer to 
the white dwarf its angular momentum must be reduced. 
However a question about the mechanism which could be 
responsible for such a reduction in the case of AE Aqr re- 
mains open. Indeed, due to a lack of a disk the canonical 
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viscosity model of the angular momentum transport is not 
applicable. On the other hand, the interaction between 
the stream and the magnetic field at a distance tq tends to 
increase the angular momentum of the material since the 
velocity of field lines, Qr, significantly exceeds the velocity 
of the material, which is limited to the free-fall velocity, 
Vs = ^IGM^^jr. 

A quantitative analysis of the mass-transfer process in 
AE Aqr has been first presented by Wynn et al. (1997). As 
they have shown, the Ha Doppler tomogram of the system 
can be reproduced assuming that the stream is inhomoge- 
neous (a sequence of large diamagnetic blobs) and interacts 
with the magnetic field of the white dwarf via a drag term. 
The efficiency of this interaction is oc r^", where n > 2. 
That is why, the strongest interaction between the blobs 
and the magnetic field occurs at their closest approach to 
the white dwarf. The initial radius of the blobs at rmin is 

Zb^lO'^Q}^' cm, (3) 

where Q\% is the efi'ective cross-section of the mass transfer 
throat at the LI point expressed in units of 10^*' cm^. The 
initial sound speed and number density of the material 
entrained in the blobs are, respectively 

Cs ^ 10*^ tI''^ cms-\ (4) 

and 

nb 2i 10^^ 9jti7 M^IJ'^ r-|/'/9 2 ^^-3^ 

Here = Zb/lO^cm, rio = rmin/lO^" cm, and is the 
mean temperature of the flowing material expressed in 
units of 10^ K (see Beskrovnaya et al. 1996; Eracleous & 
Home 1996). 

Let us now consider a scenario in which the hot fiow 
forms at Tmin as a result of strong interaction between the 
blobs and the magnetic field (independently of the mech- 
anism of this interaction). This implies a heating of the 
blobs (or their surface layers) and their expansion up to 
a size of i^. As the gas expands its density decreases by 
a factor of t\l(l\^s), where As is the thickness of the 
heated layer. This indicates that the initial number den- 
sity of material in the layer should satisfy the following 
condition 

no > nx4/(«bAs). (6) 

Otherwise, the density of the hot fiow would be signifi- 
cantly smaller than that evaluated by Itoh et al. (2005). 
Combining Eqs. (5) and (6), and setting no = rib one finds 

As > 8 X 10« OTr/ Mll^ l\, (^^t^|^) cm. (7) 

Thus, for the considered scenario to realize the entire blobs 
at a distance rmin should be heated to a temperature of 

5keV. 

However, as soon as the temperature of the blob reaches 
5 keV, its X-ray luminosity proves to be 

L,,,i ~ 6 X 10^2 t]':^ 1% ergs-i, (8) 

i.e. a factor of 60 larger that the X-ray luminosity of 
the system (~ lO^^ergs"^, see e.g. Clayton & Osborne 
1995; Choi, Dotani, & Agrawal 1999). Here riM = 
nb/lOi'^ cm-3. 

Furthermore, it remains unclear how the temperature 
of the blob could reach the value of 5 keV. Indeed, for the 
heating of the blobs to occur the characteristic time of the 



heating process should be smaller than the cooling time. 
But the bremsstrahlung cooling time of the blobs, 

ibr ^ 15 t]I^ nil s, (9) 
is significantly smaller than the free-fall time at rmin 
(iff ~ 70 t^q"^ s) and correspondingly, the charac- 
teristic time of the drag interaction, which according to 
Wynn et al. (1997) is limited to tdrag > iff- The value of 
ibr is even smaller than the turbulent diffusion time eval- 
uated by Meintjcs & Venter (2005). Thus, the assumption 
about a significant heating of the entire blobs at the ra- 
dius rmin within any currently considered mechanisms of 
interaction between the blobs and the magnetic field of 
the white dwarf is not valid. This means, that formation 
of the hot flow with the munbcr density ^ 10^^ cm^'^, lin- 
ear scale > 10^° cm, and temperature ~ 5 keV inside the 
Roche lobe of the white dwarf is impossible within the 
currently adopted model of the mass-transfer. Instead, 
the source associated with the hot surface layers of the 
blobs within this model is expected to be dense, ~ rib, 

and compact. As < 3 x 10^ nf/ T^ ^^"^ cm. 

3. DISCUSSION 

The problem posed in the previous section suggests that 
either our current view on the mass-transfer picture in 
AE Aqr is incomplete or a significant part of detected X- 
rays is generated in a region situated beyond the Roche 
lobes of the system components, or both. Here I briefly 
address these two possibilities. 

The situation can be partly improved if one invokes an 
assumption that the mass transfer through the Roche lobe 
of the white dwarf operates via more than one component. 
In particular, a formation of the X-ray source with the 
required parameters inside the Roche lobe of the white 
dwarf could be expected if the blobs were surrounded by a 
medium of a linear size ~ 10^° cm and mean number den- 
sity ~ 10^^ cm~^. The required rate of mass transfer by 
this additional flow is ~ 5 x lO^'"" g s^^, i.e. significantly (by 
more than an order of magnitude) smaller than the rate of 
mass transferred by the blobs. The heating of the fiow, ac- 
cording to Eq. (9), can be expected even if its interaction 
with the magnetic field is gentle. In particular, the cooling 
time of this fiow would be larger than the time of heating 
governed by the drag interaction (see Eq. 1 in Wynn et al. 
1997). The temperature of the flow, in the flrst approxi- 
mation, can be limited by equating the thermal pressure 
of the material, {\/2)n^m-pC^, with the energy density of 
the external magnetic field, B'^ /Sn: 

T<6.8.1ifBl {j^^Y'k. (10) 

Here B2 is the strength of the external magnetic field ex- 
pressed in units of 10^ G. 

It should be noted that both the density and size of the 
additional flow component signiflcantly exceed the corre- 
sponding parameters of the interblob material considered 
by Wynn et al. (1997). The hot flow under these conditions 
screens the blobs from the magnetospheric fleld making 
the drag interaction between them ineffective. However, 
the suppression of the drag-driven propeller action does 
not occur at the closest approach of the ffow to the white 
dwarf since the blobs and the hot ffow at this distance 
follow different trajectories. Indeed, being more dense the 
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blobs would significantly deeper penetrate into the magne- 
tosphere of the white dwarf (for a discussion see e.g. Welsh 
et al. 1998; Ikhsanov et al. 2004). This allows an ample 
room for a direct interaction between the blobs and the 
magnetic field of the white dwarf at the distance tq, i.e. 
just at the point where the efficiency of the drag interac- 
tion reaches its maximum value. 

At the same time, the origin of this additional compo- 
nent within the currently adopted model of the system re- 
mains uncertain. An incorporation of the additional com- 
ponent into the mass transfer scheme implies that either 
its temperature at the LI point exceeds 2 x 10^ K (in this 
case a formation of the geometrically thick flow at Tmin 
can be treated in terms of its thermal expansion), or the 
effective cross-section of the mass-transfer throat at the 
LI point is significantly larger than its canonical value. A 
reason why these possibilities cannot be simply discarded 
is that the normal component of AE Aqr is partly (by al- 
most a half of its radius) situated inside the light cylinder 
of the white dwarf. This situation is unique for all cur- 
rently known close binaries and the mass-transfer process 
which could be realized under these conditions is poorly 
understood so far. 

Another problem, which has to be addressed within the 
above presented scenario, is a lack of any significant or- 
bital Doppler modulation of the N and O Lya lines (Itoh 
et al. 2005). The flow effectively contributes to the X-ray 
flux of the system only on a time scale of its expansion. 



texp ^ 160 ho S, (11) 

and therefore, would appear as a local source. Indeed, 
the temperature and density of an adiabatically expand- 
ing fully ionized gas decrease as T cx Z^^, and n oc l^^ and 
the luminosity of the flow decreases on a time scale of iexp 
by more than an order of magnitude. This argues strongly 
against an assumption that the hot gas fills the Roche lobe 
of the white dwarf and spreads around the system (as it 
has been suggested by Itoh et al. 2005). Under these con- 
ditions, however, a modulation of the flow parameters with 
the system orbital motion would be expected and the lack 
of success in searching for this modulation can hardly be 
interpreted in a simple way. 

A modification of the currently adopted mass transfer 
model is not required if the source of X-ray emission has a 
multi-component nature. In particular, one can envisage a 
situation in which a part of the detected X-rays are gener- 
ated outside the system. The only available energy source 
in this case is the spin-down power of the white dwarf (the 
accretion power in this case can obviously be excluded). 
Can this energy be converted into a 5keV emission in a 
region situated outside the system? 

The answer to this question depends on the ratio of the 
spin-down power channeled into the mass outfiow and par- 
ticle acceleration respectively, which is a matter of serious 
debate at present. It is rather negative if the spin-down 
power is converted mainly into the kinetic energy of the 
outflowing material. Indeed, the velocity dispersion of the 
blobs leaving the system is limited to < 300kms~^ 
(Wynn et al. 1997; Welsh et al. 1998; Ikhsanov et al. 2004). 
This is a factor of 4 smaller than both the thermal veloc- 
ity of a plasma heated up to the temperature 5 keV and 
the velocity evaluated from the analysis of the width of 
emission lines (see Introduction). 



The answer might be, however, positive if the spin- 
down power were released in a form of accelerated parti- 
cles or /and MHD waves. As recently shown by Antonicci 
& Gomez de Castro (2005), the temperature of X-ray pho- 
tons emitted by dense (rio ^ 10^^ cm^^) blobs illuminated 
by a shower of high-energy electrons (0.03-1 Me V) ranges 
within an interval 1-20 keV. The spectrum of this radia- 
tion depends on the angle between the line of sight and the 
velocity vector of electrons. In particular, the spectrum of 
blobs situated between an observer and the source of accel- 
erated particles differs from the spectrum emitted by blobs 
situated behind the particle source. Therefore, an acceler- 
ator of particles surrounded by the blobs would appear as 
a multi-temperature source of 1-20 keV emission. 

Let us check whether such a situation can be realized in 
AE Aqr. The radius of dense blobs ejected from the sys- 
tem by the propeller action of the white dwarf increases 
due to their thermal expansion and reaches the value of h 
on a time scale 



10^ ho ce ' s. 



(12) 



Here the speed of sound, cg = Cg/lO^ cms~^, is normalized 
by taking into account that the heating of the dense blobs 
inside the Roche lobe of the white dwarf is insignificant 
(see Sect. 2). On this time scale the blobs moving with an 
average velocity V^ut leave the system to a distance of 
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VontT ~ 3 X 10" ho 



out 



300kms-i 



cm 



(13) 



and their number density decreases to 
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no = nb(^]'^10" ll £^o'{ 



1014 cm-3 



) cm-3 (14) 



The number of blobs ejected by the white dwarf during 
the time r is 

^ ^ 3 grtr 
~ ATrnhrUpll ' 

Combining Eqs. (5), (12), and (15) yields 



N ~ 10^ Mo'/' Cg 1 ho /g ' 



-1/2 



10 • (16) 

Therefore, the stream of the expanded blobs will occupy 
an area of 



A ~ welN ~ 3 X 10' 



23 



4 



1010 cm 



TV 

103 



cm 



(17) 



which constitutes ~ 0.3 of the total area of a sphere with 

a radius of tq. This indicates that 30% of the energy re- 
leased in the system will be transferred through the stream 
of the expanded blobs. 

If now we assume that the spin-down power releases 
predominantly in a form of accelerated particles (see e.g 
Ikhsanov 1998; Meintjes & de Jager 2000) one finds the 
total flux of the energy transferred through the area occu- 
pied by the stream of the expanded blobs as (see Eqs. 5, 
and 12-17) 



L, =L. 



sd 



47rro 



~ ergs-i hoMl's ce Iq^ ho'^ 



(18) 



300 kms- 



L. 



sd 



6 X 10^3 ergs~i 



4 



Hence, a conversion of only 2% of this energy into X-rays 
would be enough for the himinosity of the corresponding 
source to be comparable to the X-ray luminosity observed 
from the system. 

The above estimates suggest that an interpretation of 
the X-ray emission of AE Aqr in terms of expanded blobs 
illuminated by a flux of accelerated particles might be 
promising. A further analysis of this scenario, however, 
requires information about the spectrum and geometry of 
the wind of particles ejected by the white dwarf, which 
is currently not available. Furthermore, a question about 
the region of the wind formation also remains open. If 
this region (as in the case of spin-powered pulsars, see e.g. 
Beskin, Gurevich & Istomin 1993) is situated at the light 



cylinder, the surface of the normal companion would also 
be partly affected by the wind. However, a contribution 
of this additional source into the X-ray flux can unlikely 
be significant since the area occupied by the normal com- 
panion at the radius of the light cylinder is a factor of 20 
smaller than A (see Eq. 17). 
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